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The 



Semiclassical theories of electron and spin transport in metallic magnetic structures are reviewed 
with emphasis on the role of disorder and electronic band structures in the current perpendicular 
to the interface plane (CPP) transport configuration. 
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I. INTRODUCTION 

Electron transport in layered magnetic systems can be directed parallel or perpendicular to the interfaces. The 
physics of transport in the latter, so-called current perpendicular to the plane (CPP) configuration, studied first by 
the Michigan State University Group with superconducting contacts 0] and by Gijs c.s. in microstructurcd pillars 
at arbitary temperatures Q , has been reviewed quite recently |^ . The topic remains to attract the interest of the 
community and new insights have been obtained on issues as first-principles calculations of transport in disordered 
multilayers, phase coherence effects, transport with non-coUinear magnetizations and spin-torques induced by applied 
currents. The present manuscript briefly reviews and compares these novel developments. Related topics of interest 
are superconductor- ferromagnet hybrids |^,^, spin- injection into carbon nanotubes |^,^ and semiconductors |^ or 
many-terminal devices Q . 

Most studies have been carried out on multilayers consisting of many bilayer periods. In earlier studies attention 
was focussed on coUinear magnetization profiles, i.e. that with antiparallel or parallel magnetization vectors for 
neighbouring layers. These experiments are well fitted by the so-called 2-channel series resistor model (2CSRM) 
which is discussed in the light of recent developments in Chapter 2. In Chapter 3 novel theoretical approaches are 
reviewed which describe non-collinear configurations. 



II. COLLINEAR MAGNETIZATION 



o 



The experiments of perpendicular transport in magnetic multilayers are well described by the 2CSRM, i.e. an 
equivalent electric circuit of two spin channels in parallel, in which the resistance of each channel is the sum of the 
bulk and interface resistances. When all magnetizations are parallel the total resistance Rt of a ferromagnetic/normal 
metal (F/N) multilayer of M double layers reads 
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where A is the cross section of the sample, p^^^ is the resistivity of the bulk normal metal, pf^ is the resistivity of 
the ferromagnet for spin direction s, diq and dp are the layer thicknesses and R'^^^ are the spin-dependent interface 
resistances. The five parameters p'-^-', pi^\ and R^^^ can be determined accurately by fitting experiments on 
numerous samples with different layer thicknesses and for antiparallel as well as parallel magnetic configurations [0. 
It was soon realized that the microscopic basis was to be found in semiclassical arguments [^,0. The most flexible 
theoretical framework turned out to be the linearized Boltzmann equation in the relaxation time approximation by 
Valet and Pert [Q, which also included spin- flip processes. In this model the distribution function or local chemical 
potential fs for spin s in the F ot N bulk materials is governed by the one-dimensional diffusion equation: 
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dx^ ~ 2 ■ ^ > 

where as is the conductivity and spin-flip processes are included in terms of the spin-diffusion length Isj . Interfaces 
were treated as thin regions with a different (low) mobility. The current for spin s then reads: 

Js{x)=(Ts , (3) 

When Igf ^ d the results are identical to the 2CSRM. However, this treatment is incomplete in that the discontinuities 
in the electronic structure at heterointerfaces are disregarded, which are essential for the electron transport properties 
since they scatter electrons even in clean samples This can be seen most easily for an isolated interface in a 

constriction according to the Landauer-Biittiker formula : 

1 _ 
R^/B ~ h 
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where t^^>,T^^i are the transmission coefficients and probabilities of states fj, — fcyi^s at the Fermi energy with 
(tranverse) wave vector fc|| parallel to the interface and band index v. For a ballistic interface fcy and s are conserved 
during scattering. This expression can be compared with that for a homogeneous point contact of the materials A 
and B with (Sharvin) conductances G'^ > . It is clear that any mismatch in the electronic structure will reduce 
the conductance G^^^ < G^. 



A. Interface resistance 



Of interest is the microscopic explanation of the interface resistance parameter R^^^ in Eq. (|l]) and its relation with 
Eq. (^). The success of the 2CSRM provides the guidance for a quantitative understanding of the experiments. The 
parameters which fit so many experiments turn out to be universal for a given material combination. The absence of a 
measurable dependence on the geometrical parameters and dp (when everything is kept the same) is a strong indi- 
cation that quantum interference terms are negligibly small. For the large thicknesses which have predominantly been 
studied experimentally, this is not surprising. But also for significantly exchange-coupled multilayers, quantum well 
states should not significantly affect transport since disorder causes a large semiclassical background. A semiclassical 
model therefore should be appropriate for all but the cleanest samples with very thin layers. A simple yet efficient 
first principles procedure to incorporate the interfaces comes down to chopping the sample into slices of bulk layers 
and interfaces which scatter electrons, separated by fictitious non-scattering regions [^. The bulk layer scattering 
can be treated in terms of transmission and reflection matrices, modelled in two simple limits; the ballistic limit, in 
which no scattering occurs during the transmission through the bulk and the diffuse limit in which the scattering is 
isotropic, which means that the transmitted electrons do not retain any memory of the incident wave vectors. The 
transmission probabilities for the combined system then follows by the semiclassical concatenation of transmission 
and reflection probability matrices of the resistive elements in series [^6|,Q. Indeed, the 2CSRM is recovered with 
interface resistances 
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where we have corrected for the spurious geometrical Sharvin resistances 



(5) 



L8|. For ballistic bulk layer transmission 
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which is in fact the "old" Landauer formula. Both Eqs. (^ and ^ do not contain any parameters and provide a 
first-principles prescription for the semiclassical interface resistances. They hold not only in the metallic case but 
should be also valid when the interface is in the tunneling regime. Stiles and Pcnn |l^] confirmed that the expression 
(^ agrees very well with the solutions of the Boltzmann equation. 

Levy et al. pO| recently found strong effects of disorder in the leads on transport through tunnel junctions. Our 
discussion implies that these effects should be absent in a semiclassical description, which is at odds with the conclu- 
sions of Levy et al. Interference effects between impurities and resonance states at the tunnel junction might provide 
an explanation [ pT| . 
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B. First-principles calculations 



The interface resistance has been calculated by first principles for specular interfaces in Refs. |15| , |l9| , p3| for ballistic 
magnetic domain walls in Ref. and diffuse interfaces in Ref. In the table we summarize our results for 

specular and rough interfaces, the latter modelled as a 50%/50% bilayer alloy, and compare them with available 
experiments. Remarkable is the different behaviour of the Co|Cu as compared with the Fe|Cr interface with respect 
to the interface roughness: the spin contrast of the Co|Cu interface transparency is weakly enhanced by disorder, 
but strongly reduced in Fe|Cr. The message for experimentalists is that in Fe|Cr it should pay off to optimize the 
epitaxial growth parameters. 

The Boltzmann equation has been solved numerically by Butler et al. [ p6[ for Co|Cu|Co perpendicular spin valves 
(see also [p7|). The constant relaxation time approximation for the bulk materials was used and distribution functions 
were matched via the transmission and reflection coefficients for specular interfaces. They found results for the 
interface resistance slightly different from that of . Butler et al. find a dependence of the interface resistance on 

small Cu thicknesses, which appears to contradict the universality of the parameters of the 2CSRM. These corrections 
are not the small quantum size effects found by Xia et al. ||2^ and by Tsymbal et al. [^,|5| , but are purely classical 
effects due to evanescent terms in the distribution functions. These can be interpreted as the corrections to the 
assumption of complete isotropy by Schep et al. [|l5| . The isotropy conditions is expected to be much better fulfilled 
when interface disorder is taken into account. It is clear that more, also experimental, work is needed to understand 
transport for very thin layer thicknesses in which the incomplete randomization and the appearance of quantum 
corrections will cause deviations from the two-channel resistor model. 



III. NON-COLLINEAR MAGNETIZATION 



Electron transport in devices with non-coUinear magnetizations have come into focus by the recent interest in the 
torque exerted on the magnetization by a spin-polarized injected current . The theories for collinear magnetization 
mentioned above have been extended to the case in which the magnetizations are not collinear, i.e. not parallel or 
antiparallel, namely via a magnetoelectronic circuit theory |29[| , and for CPP spin- valve structures, random matrix 
theory of transport pof as well as a direct solution of the diffusion equation in the presence of an external magnetic 
field |3^|. Recently also exchange-biased CPP spin- valves are under scrutiny Bo]. We find that these approaches are 



equivalent and reduce to previous theories in the limit of collinear magnetization. An interesting point is that ^9 30 
do not start from the outset with a semiclassical approximation, but derive it from an isotropy assumption, thus put 
previous more or less ad hoc approaches on firm theoretical foundations. 



A. Circuit theory 



Transport in hybrid metallic systems can be described by a generalization of Kirchhoff 's theory of electronic circuits 
when parts of the system are not phase-coherently coupled. This approach has been pioneered in Ref. [ |3^ for electronic 
networks with superconducting elements. It has recently been adopted also for magnetoelectronic circuits like 
the Johnson spin transistor |Q| or the 4-terminal mesoscopic spin valve of Jedema et al. . For a different approach 
to the many terminal magnetolectronic circuits see |^,^. The circuit theory can be derived from a given Stoner 
Hamiltonian in terms of the Keldysh non-equilibrium Green function formalism in spin space . The basic physics 
is provided by splitting up the system into reservoirs, resistors and nodes, where the latter can be real or fictitious 
(as discussed above). In order to arrive at a useful formalism, an isotropy assumption has to be introduced for the 
nodes, namely that the electron distributions in the nodes are isotropic, which implies sufficient disorder (or chaotic 
scattering) to allow for configurational averaging. It does not require any inelastic or dephasing scattering mechanism, 
although, when happening in the nodes, they will not hurt cither. Because the spin- accumulation is not necessarily 
parallel to the spin-quantization axis, at each node the electron distribution at a given energy e can be denoted as /(e), 
where the hat denotes a 2 x 2 matrix in spin-space. The external reservoirs are assumed to be in local equilibrium 
so that the distribution matrix is diagonal in spin-space and attains its local equilibrium value / = l/(e, ^J■a)l 1 is the 
unit matrix, f(e,fia) is the Fermi-Dirac distribution function and fia is the local chemical potential in reservoir a. 
The direction of the magnetization of the ferromagnetic nodes is denoted by the unit vector rxia . The current through 
each contact can be calculated as a function of the distribution matrices on the adjacent nodes 2x2 conductance 
tensors. The current matrix (for an F|N junctions) reads 

ei = GT^T (^fP _ /A'j + G^ii^ (/^ - /^) - G^^u^f^u^ - (G^^yu^f^'u^ , (7) 
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in terms spin-irotation matrices = (I + sa ■ mj /2, distribution matrices y j'^ on ferromagnetic and normal 
node, the spin-dependent conductances and (which in planar junctions should be corrected as discussed in 
Section II) 



and the mixing conductance 



G" = — 
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(8) 



where r"™, t"™ are the reflection and transmission coefficients, M the number of modes in the absence of reflections. 
Spin-flips in the contacts have been disregarded. The spin-current conservation law 



^ dU\ _ ITr/^ - 2/^ 



dt 
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(9) 



allows computation of the circuit properties as a function of the applied voltages, where lap denotes the 2x2 current 
in spin-space from node (or reservoir) a to node (or reservoir) /3 and the term on the right hand side describes spin- 
relaxation in the (normal) node. The right hand side of Eq. m) can be set to zero when the spin-current in the node 
is conserved, i.e. when an electron resides on the node sufficiently shorter than the spin-flip relaxation time r^j. 

Some insight can be gained by re-writing the current and the distribution function in the form of a scalar particle 
and a vectorial spin contribution, I — {Io + (t-1s)/2, — f.^ +(t- sAf^ and j — fp +(t ■ mAf^ . The spin-current 
through an F|N interface can then be expanded into different vector components as: 



I. - m[(GT - G^){f^ ~ /p^) + (GT + G^)Af + (G^ + G^ 



2ReG^^)s-mA/^] 
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The vector spin current component perpendicular to the magnetization direction equals the spin-torque exerted by 
the polarized current on the ferromagnet |2q,pG]. 



B. Random matrix theory 



Waintal et al. \ 30 1 have extended the random matrix theory of transport |^ to include non-coUinear magnetizations 
in F|N|F junctions. This paper is focussed on the spin torque, but the physics is essentially the same as for transport 
in a disordered CPP spin valve with arbitrary magnetization configuration as discussed above. The objects which are 
averaged are the scattering matrices of the bulk layer of the normal metal assuming that all members of the ensemble 
fulfil the symmetry requirements of the problem and are equally probable. The theory is more intricate than in the 
case of superconducting S|N|S junctions because the averaging has to be carried out over the eigenvalues and 
eigenvectors. Analytical results can be obtained for the leading term of an expansion into l/N , i.e. the inverse of the 
number of transport channels. 

It is easily seen that the analytical relations obtained by Waintal et al. |Q for halfmetallic ferromagnetic elements 
agree with those from the circuit theory for the symmetric two-terminal device. The general case is less obvious, but 
by somewhat tedious manipulations a complete equivalence of the final equations for both theories can be proven 
psf . We may conclude from Waintal's results that what seemed to be an assumption in the procedure of Schep et 
al. , viz. the semiclassical concatenation, can be derived from the isotropy assumption. Although not yet worked 
out, Waintal et al. 's approach can be generalized to include quantum corrections, which become important for a small 
number of channels, as well as many terminal configurations. It is not clear how spin-flip-relaxation processes can be 
incorporated, which is quite straight-forward for the circuit theory. 



C. Diffusion equation for non-collinear transport 

In some cases the theories above are not sufficient and the spatially dependent ditribution has to be evaluated. The 
spin-polarized electron distribution is characterized by a 2 x 2 matrix in spin space of the form: 
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(11) 



When the size of the system L is larger than the spin diffusion length Isf, f^{x) depends on the position. We have 
studied transport through an F|N|F device under the condition If Isf, where If = u_f (f /r + is the mean 

free path, vp is the Fermi velocity, r the spin-conserving scattering time and r^/ the spin- flip scattering time [ pf| . 
Under the condition / f <C hf — \/^fI fTsf /3, we obtain the generalized diffusion equation in the normal metal 
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(12) 



Its solution, with boundary conditions at the interface governed by the conductance matrix, describes e.g. the 
precession of the spin-accumulation in an applied magnetic field and lead to a physical interpretation of the imaginary 
part of the mixing conductance [pT| . 



D. CPP spin valve 

The different approaches described above lead to an analytical expression for the total conductance of CPP spin 
valves as a function of the angle between the magnetizations of the different ferromagnets 9, when l^f ^ L, at zero 
magnetic field {B = 0) and for symmetric contacts: 
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where G ^ + , P = (G^ + G^) /G, 



_2G^ hi" _ 4|GTi|2 Q ^ 2|GTi|2 

" ~G~'' ~ G2 2ReGn ~ GReGTi ^^^^ 
The angular magnetoconductance reads: 

G^(g) -G^(0) _ tan2g/2 _ [ sin" 0/2 r, = Rer,^l 

G^in) - G-(0) - tan^ 9/2 + ^^' ^''^ 

Balents and Egger arrive at the same result in their study of spin-injection into carbon nanotubes in the non- 
interacting limit. The electron-electron interaction is found to enhance the mixing conductance. 



IV. DISCUSSION 



Our understanding of the transport properties of the CPP multilayers is semi-quantitative for the parallel aligned 
and the as grown "virgin" samples, in which to a good approximation neighbouring magnetization vectors are antipar- 
allel. The basis of this understanding is (I) knowledge of the magnetization configuration and (2) the 2CSRM. Bozec 
et al. [Q claimed to have found evidence for a breakdown of the 2CSRM. However, a recent study with intentionally 
alloyed bulk layers comes to different conclusions, i.e. that the 2CSRM should remain unchallenged for coUinear 
magnetic structures, be they "type I" or "type H", interleaved or separated [^ . 

The situation of the magnetic- field cycled "deflowered" samples is more difficult. The experiments of Bozec et al. 
could be explained by a "spin-memory" effect caused by spin-flip at the interfaces which can be incorporated 
into the two-channel series resistor model [ p9[ . This picture requires that the magnetization at all intermediate fields 
is random but essentially coUinear (except possibly at interfaces). An alternative explanation is Wiser's hypothesis 
that the angle between the magnetizations of different layers is rotated during magnetization reversal |^^. The 
transport properties based on this hypothesis can be computed in principle by the generalizations of the 2CSRM 
to non-coUinear transport discussed above, which we may call "matrix series resistor model" . At the moment the 
magnetization distribution is not known sufficiently well, but it seems likely that non-coUinearity and randomness 
both play a role. Exchange-biased CPP spin valves appear to be better suited to test the new theories than multilayers 
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system 


roughness 


Rmaj(fnni2) 


Rmin(fnm^) 


Au/Ag(lll) 


clean 


0.094 


0.094 


Au/Ag(lll) 


2 layers 50-50 alloy 


0.118 


0.118 


Au/Ag(lll) 


exp. i 


0.100 ± 0.008 


0.100 ± 0.008 


Co/Cu(100) 


clean 


0.33 


1.79 


Cohcp/Cu(lll) 


clean 


0.60 


2.24 


Co/Cu(lll) 


clean 


0.39 


1.46 


Co/Cu(lll) 


2 layers 50-50 alloy 


0.41 


1.82±0.03 


Co/Cu(lll) 


exp. § 


0.26±0.06 


1.84±0.14 


Fe/Cr(100) 


clean 


2.82 


0.50 


Fe/Cr(100) 


2 layers 50-50 alloy 


0.99 


0.50 
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